Quasiparticle States around a Nonmagnetic Impurity in D-Density-Wave State of 

High-T c Cuprates 
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Recently Chakravarty et al. proposed an ordered d-density wave (DDW) state as an explanation of 
the pseudogap phase in underdoped high-temperature cuprates. We study the competition between 
the DDW and superconducting ordering based on an effective mean-field Hamiltonian. We are 
mainly concerned with the effect of the DDW ordering on the electronic state around a single 
nonmagnetic impurity. We find that a single subgap resonance peak appears in the local density of 
state around the impurity. In the unitary limit, the position of this resonance peak is always located 
at E r — — fi with respect to the Fermi energy. This result is dramatically different from the case of 
the pure superconducting state for which the impurity resonant energy is approximately pinned at 
the Fermi level. This can be used to probe the existence of the DDW ordering in cuprates. 



PACS numbers: 74.25Jb, 74.20.-z, 73.20.Hb 

Recently a d-density wave (DDW) state was pro- 
posed to model the pseudogap in the underdoped high-T c 
cuprates (|||| . The key feature of the DDW state is stag- 
gered orbital magnetic moments (i.e., staggered currents) 
which break parity and time-reversal symmetry as well 
as the translation invariance by one lattice constant and 
rotation by tt/2. This proposal of a new order parameter 
differs fundamentally from other theoretical ideas |^-||] 
for which the staggered current is a fluctuating rather 
than a static quantity. Whether such a DDW state re- 
ally exists in the phase diagram of high-Tc cuprates is an 
interesting theoretical question Q. More recently, the 
SU(2) § or U(l) § mean-field theory of the t-J model 
predicted that the DDW order could exist in a vortex 
core. 

An interesting question is the effect of a single non- 
magnetic impurity on the DDW state. The single im- 
purity problem in the superconducting state has been 
intensively studied both theoretically [pHl2|] and exper- 
imentally |l3|-|l5[] ■ However, the same problem has not 
been well studied in the pseudogap (PG) state, which 
has been indicated by STM |1J and by ARPES exper- 
iments jl7j. Currently, the origin of the PG state re- 
mains unclear. Experimentally it was also argued |is[] 
that in the underdoped cuprates the superconducting 
and pseudo gap coexist unlike in the pre-formed pair 
model The electronic states around a single im- 

purity may serve as a local probe to distinguish between 
these scenarios. Within the T-matrix approximation, it 
was argued Q that, as long as the band density of states 
is depleted at the Fermi energy, the existence of impurity 
resonant states is robust regardless of the microscopic 
origin of the PG state. Note that, without invoking any 
particular model, the band density of states was assumed 
in Ref. [^0| to be symmetric with respect to the Fermi en- 
ergy. This choice corresponds to the chemical potential 
/i being zero in a two-dimensional tight-binding model 
with nearest-neighbor hopping. In addition, because of 



the lack of a specific model, the calculation of the local 
density of states (LDOS) around an impurity goes be- 
yond the scope of that work. The purpose of this paper 
is to address this question within the DDW state, for 
which the LDOS around the impurity is calculated. Our 
main result is that, in the strong scattering limit, the 
position of the induced resonant state has a one-to-one 
correspondence to the chemical potential, which is dra- 
matically different from the case of the superconducting 
state (and also of the pre-formed pair scenario), where 
the resonant energy is almost pinned at the Fermi en- 
ergy. This prediction can be readily used by the STM 
experiments to detect for the DDW state in the under- 
doped cuprates. 

We start with a phenomenological model including the 
intersite charge density wave order in a d-wave super- 
conductor. One could also envision intersite spin density 
wave instead of the charge density order; however, as long 
as the density order is to be interpreted as a pseudogap 
the charge density order is more reasonable experimen- 
tally M. The effective Hamiltonian is written as: 

H = ]p-i + (-l)Wy]4c jV + J2(Ui - M )cU CT 

ij.cr i.cr 

^(Ayc^+A^cit). (1) 



Here c ia creates an electron with spin a at site i. The 
summation is over the nearest neighbors sites, t is the 
hopping integral and fi is the chemical potential. U\ is 
introduced to model the potential scattering, if any, from 
impurities or defects. Wi- } and Ay are, respectively, the 
DDW and DSC order parameters, which are determined 
by the self-consistent conditions: 



^ = (-1)^(4^-4^ 



(2) 



and 
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(3) 



where Vd dw and Vdsc? are the interaction for the DDW 
and DSC channels. (— l) 1 in the ansatz Eq. (|J) ensures 
that nesting takes place at half filling, and that Wi- } is 
pure-imaginary. A. Ghosal et. al. |2l]] have considered a 
Fock shift, which is similar with W/y but does not have 
the site-dependent factor (— 1)". A possible model for 
the Hamiltonian |l| is the t-J model. Applying the mean- 
field approximation and the ansatz for the density wave 
to the exchange interaction, one can obtain density and 
superconducting order. Introduction of an intersite in- 
teraction in addition to the exchange interaction gives 
different values of Vdsc and Vddw- 

The mean-field Hamiltonian ( P can be diagonalized 
by solving the resulting Bogoliubov-de Gennes equations 
self-consistently 



vV 



(4) 



where the single particle Hamiltonian TL-^ 



-td 



(— l) 1 Wij + (Ui — fJ-)Sij . and the self-consistency now reads 



J2 Im[<up* + vl'vf*} tanh 



\2k B T J 
(5) 



and 



Competition between DDW and DSC orderings. The 
phase diagram with the DDW and DSC orderings in 
Ref. [pj was constructed from experiments. Therefore, 
it is important to self-consistently study the phase di- 
agram of a clean system with DDW and DSC channel 
interactions. Although the DDW state breaks the trans- 
lational invariance by one lattice constant a, it is in- 
variant under the translation with spacing y/2a along 
the diagonal directions of the square lattice. Hereafter 
we will measure the length in units of a and the en- 
ergy in unit of t. The site-dependent factor (— l) 1 en- 
ables us to consider the square lattice as a bipartite lat- 
tice. Choosing the basis cell with the size 
we solve the BdG equation ( |J) and find the eigenval- 
ues E k = ±£i )2 (k) with £ 1)2 (k) = [(& T M) 2 + A£] 1/2 , 
where e k = — 2(cosfc K + cosk y ), Wk — 2iWd{cosk x — 
cosky), £ k = V e k + \ w k\ 2 and A k = 2A (J (cosfc a; - 
cos k y ). The eigcnfunctions corresponding to (k) are, 
(u^e^ DDW ,±Uk,Vke~^ DSC ,±vi Q e~ i ^ DSC ~^ DDW ' ) ) Transpose . 

Here u k = ± (l ± and «M(lT , and 

4>ddw = tan _1 (— iWk/ek) and 4>dsc are the phase as- 
sociated with the DDW and DSC orders. Substitution 



of these eigenfunctions into the self-consistent equations 
for DDW and DSC order parameter leads to: 



W d 



Vddw W d (co$k x - cosfc y ) 2 £ k + (-l) Q /x 



E 



x tanh ■ 



2k B T 



Vdsc A d (cosfc a; - cosky) 2 



k - 



E a 



2k B T ' 



(7) 



(8) 



where a — 1,2 and Nl is the number of the basis cells. 
Notice that at the fixed doping, the chemical potential 
itself is determined from 5 = 1 — n where the filling factor 
n =WI Ek,J/(£aK + [1 - f(E a )H}. The solutions 
to the above equations have shown that the critical tem- 
perature T DDW for the "bare" transition to the DDW 
state decreased more quickly than T DSC for the DSC 
state. Therefore, only when the DDW channel interac- 
tion is larger than the DSC one, can there exist separated 
regions where either the DDW or DSC state is dominant. 
As a model calculation, we display in Fig. [I] the temper- 
ature versus doping phase diagram for Vddw = 1.6 and 
Vdsc — 1-4. In the region of doping close to zero, the 
DDW state is dominant while in the region of S away from 
zero, the DSC state is dominant. In addition, we have a 
finite area of coexisting DDW+DSC phase. The size of 
the pure DDW region could be increased with a differ- 
ent choice of Vddw and Vdsc but this is not expected 
to lead to important quantitative change. Also half fill- 
ing in our model is not an insulating anti-feromagnetic 
phase because we have not included any of this physics 
in our model so that the pseudogap state extends down 
to S = 0. 

This competition can be understood from a phe- 
nomenological Ginzburg-Landau (GL) theory. The GL 
free energy density, in terms of both the DDW and DSC 
order parameters, can be constructed from a symmetry 
analysis: 



/ = a D sc\&d\ 2 + a DD w\W d \ 2 + ft | A d 



P2\W d \ 2 



+P 3 \A d \ 2 \W d \ 2 + M&fW* + A 2 Wf) , (9) 



where we assume oidsc 



1 dsc 



(T 



Tdsc) 



and 



c*ddw = a DDw(T — T DDW ), and f3's are all positive. As 
an example, in the region where T DDW > T DSC , the sec- 
ond phase transition temperature for the DSC ordering is 
rcnormalized by the pre-existing DDW order parameter: 



1 DSC = J-DSC ~ 777, } 

2ha, 



DSC 



DSc/ a DDW 



(ft - 20 4 ) 



(10) 



Therefore, even if the "bare" critical temperatures for 
both ordering are very close to each other, the second- 
phase transition temperature for the appearance of DSC 
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(or DDW) order parameter can be strongly suppressed by 
the dominant DDW (or DSC) order parameter. At the 
mean- field level, we conclude that when Vddw > Vdsc-, 
the phase diagram is consistent with that proposed in 
Rcf. 0| . Here we also would like to point out that the 
U(l) mean-field theory of the t-J model gives rise to 
Vddw = 0.5J while Vdsc — J, which may explain the 
absence of the bulk DDW state in that model |(| . 

The LDOS around a single nonmagnetic impurity. The 
LDOS is given by: 

Pi (E) = -2j2M\ 2 f(E n -25) + \v?\ 2 f'{E n +E)] , 

n 

(ii) 

where a factor 2 arises from the summation over spin, and 
f'(E) = df(E)/dE with the Fermi distribution function 
f(E) = l/{e E / kBT + 1). The LDOS p v {E) is proportional 
to the local differential tunneling conductance which can 
be measured by STM experiments (22) . In our numerical 
calculation, we take the supercell size N x = N y = 32 and 
the number of supercells M = 6 x 6. Since we are most 
interested in identifying the qualitative difference in the 
electronic states around the impurity for the DDW and 
DSC states, the thermally broadening effect will not be 
considered here and the temperature is fixed at T = 0.01. 
For simplicity, we also use uniform order parameter as an 
input to diagonalize Eq. ([IJ) and the suppression of the 
order parameter near the impurity is ignored. This ap- 
proximation should be acceptable to answer the questions 
we ask here. 

We first consider the case of weak or moderately strong 
impurity scattering. Figure || gives a plot of the LDOS 
spectrum directly on the single-site impurity (a) and on a 
site one lattice constant away (b) for a pure DDW state. 
When the impurity scattering is weak (Uq — 4), the 
LDOS on both the impurity site and on its nearest neigh- 
bor has a single peak below the Fermi energy. Although 
the LDOS at the impurity site is similar to the case of a 
weak nonmagnetic impurity in a pure DSC state |]l0| , ^3| , 
the LDOS displays a single resonant peak at the nearest 
neighbor site while the DSC shows a double-peak struc- 
ture with one peak above and the other below the Fermi 
energy. When the impurity scattering becomes stronger 
(Uq = 100), the subgap resonant peak in the LDOS at 
the impurity site is shifted toward E r = with respect 
to the Fermi energy at the same time its amplitude at 
the impurity site is strongly suppressed because of the 
strong impurity scattering while at the nearest neighbor 
site it is strongly enhanced. 

In the unitary limit, the energy position E r for the 
impurity resonant state in a pure DDW state is very sen- 
sitive to the chemical potential. In Fig. ||, the LDOS 
spectrum is plotted at the nearest-neighbor site for var- 
ious values of fi. We see in Fig. |^ that the resonant 
energy position is exactly equal to E r = (p < 0), 



which is in sharply contrast to the case of a single im- 
purity in a pure DSC state, where the resonant energy 
is not sensitive to the chemical potential and roughly 
zero (i.e., very close to the Fermi surface). This differ- 
ence can be understood as follows: The energy dispersion 
for the pure DDW state £a, 2 (k) = \y/e£ + |!F k | 2 T A«l 
shows that the overall quasiparticlc band is shifted by fi 
because of nature of the DDW state. Correspondingly, 
the resonant peak which always exists at the band cen- 
ter is shifted by /i. However, due to the pairing mecha- 
nism, the quasiparticle for the pure DSC state, as given 
by Ei i2 (k) — [(|ek| T m) 2 + Ak] 1 ^ 2 : are always excitated 
w.r.t. the Fermi surface instead of the band center. Con- 
sequently, the impurity induced resonant state has the 
energy almost pinned around the Fermi surface. More- 
over, as in the case of a DSC state, we have found that 
the impurity resonant peak intensity of the LDOS in a 
pure DDW state exhibits the Friedel-like spatial oscilla- 
tion: It has local maxima on the sublattice containing 
the nearest neighbors and local minima on the sublattice 
containing the impurity site itself. 

Finally, in Fig. ||, the LDOS spectra at the nearest- 
neighbor site of a moderately strong impurity is plot- 
ted in a mixed state of DDW and DSC ordering. For a 
comparison, the LDOS spectrum at the same site is also 
displayed in a pure DSC state. In the pure DSC state, 
we have the double-peak structure with the intensity of 
the peak above the Fermi energy stronger than that be- 
low the Fermi energy. However, in the mixed state with 
both orderings present, the intensity of the peak below 
the Fermi energy is stronger than that above because the 
impurity scattering in the DDW order shifts the resonant 
peak below the Fermi energy. 

In conclusion, we have studied the competition be- 
tween the DDW and DSC ordering. The implication of 
the DDW state for the impurity resonant state has been 
discussed in detail. The qualitative difference found in 
the resonant state in the DDW ordering compared with 
that in the DSC state can be used as a smoking gun 
for the existence of the DDW state in high-T c cuprates. 
Experimentally, the existence of this state can be identi- 
fied by (i) detecting whether the local differential tunnel- 
ing conductance near a single weak or moderately strong 
nonmagnetic impurity exhibits a single rather than dou- 
ble subgap structure around the Fermi surface and (ii) 
detecting the position of the subgap resonant peak in 
the conductance near a unitary nonmagnetic impurity 
which will have a strong doping dependence. An STM 
measurement is most suitable for this test. 
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FIG. 3. The LDOS spectrum at the nearest neighbor of a 
unitary impurity in a pure DDW state for various values of the 
chemical potential /i = (red-solid line), —0.1 (green-dashed 
line), and —0.2 (blue-dash-dotted line). The other parameter 
values: W d = 0.1, U = 100, and T = 0.01. 



FIG. 4. The LDOS spectrum at the nearest neighbor of a 
weak impurity in a mixed state of DDW and DSC orderings 
(red-solid line). The other parameter values: Wd = 0.08, 
A d = 0.1, Uo = 4, fi = 0, and T = 0.01. Also shown is the 
spectrum in a pure DSC state with A d = 0.1 (green-dashed 
line). 



FIG. 1. The temperature versus doping(<5) phase diagram 
obtained from the Hamiltonian (|l|) with channel interaction 
Vddw = 1.6 and Vdsc = 1.4. 



FIG. 2. The LDOS spectrum at the impurity site (a) and 
at its nearest neighbor (b) in a pure DDW state for various 
strength of the impurity scattering Uo = 4 (red-solid line) 
and Uo — 100 (green-dashed line). Also shown the LDOS 
spectrum (blue-dash-dotted line) at a position far away from 
the impurity. The other parameter values: Wo = 0.1, /i = 0, 
and T — 0.01. The blue line in the bulk density of states for 
reference. 
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FIG.1 Zhu etal. 




FIG.2 Zhu et al. 




FIG.3 Zhu et al. 




FIG.4 Zhu et al. 



